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Plaquette lattices with each unit cell containing multiple atoms are good candidates for discon- 
nected Fermi surfaces, which are shown by Kuroki and Arita to be favorable for spin-flucutation 
mediated superconductivity from electron repulsion. Here we find an interesting example in a 
tetramer lattice where the structure within each unit cell dominates the nodal structure of the gap 
function. We trace its reason to the way in which a Cooper pair is formed across the hybridized 
molecular orbit als, where we still end up with a T c much higher than usual. 

PACS numbers: 74.70.Kn, 74.20.Mn 



The discovery of high-T c superconductors has 
kicked off renewed interests in electronic mechanisms for 
superconductivity. In the course of studies it is becom- 
ing increasingly clear that superconductivity can arise 
from repulsive interaction between electrons. The essence 
there is electrons interact by exchanging spin fluctua- 
tions, which can effectively act as an attraction in the gap 
equation, if we consider anisotropic pairing with nodes in 
the gap function. 

For the Hubbard model, this was first suggested from 
early calculations by Scalapino et al. [||. Subsequently 
a quantum Monte Carlo calculation || indeed indicated 
an enhancement of the pairing correlation with <i-wave 
symmetry in the repulsive Hubbard model. The super- 
conducting critical temperature (T c ) has been estimated 
with the fluctuation-exchange approximation (FLEX), a 
kind of renormalized random phase approximation. Q-|4|] 

One remarkable point is T c ~ O(0. Olt), esitmated 
for the repulsive Hubbard model in the two-dimensional 
(2D) square lattice, is two orders of magnitudes smaller 
than the starting electronic energy (i.e., the hopping inte- 
gral i), although this gives the right order for the curates' 
T c . Arita et al. then looked at various lattices (square, 
trianglar, fee, bec, etc) in search of the most favorable 
case with a FLEX analysis 0. The best case, as far 
as these ordinary lattices are concerned, turns out to be 
the 2D square lattice, so T c < O(0. Oli) remains. As dis- 
cussed in Ref. ||, there are good reasons why T c is so low, 
where an important one is the presence of nodes in the 
superconducting gap function greatly reduces T c : While 
the main pair-scattering, across which the gap function 
has opposite signs to make the effective interaction at- 
tractive, some of the pair scatterings around the nodes 
have negative contributions to the effective attraction by 
connecting /c-points on which the gap has the same sign. 

So a next important avenue to explore is: can we 



improve the situation by going over to multiband sys- 
tems. Kuroki and Arita have shown that this is in- 
deed the case if we have disconnected Fermi surfaces. In 
this case T c is dramatically enhanced, because the sign 
change in the gap function can avoid the Fermi pockets, 
where all the pair-scattering processes contribute posi- 
tively. p|-|Io[ This has been numerically shown to be the 
case for the triangular lattice (for spin-triplet pairing) 
p| and a squre lattice with a period-doubling H , where 
T c as estimated with FLEX is as high as O(0. It). 

To be more precise, the key ingredinents are: (a) when 
the Fermi surface is nested, the spin susceptibility x(q, w) 
has a peak, with a width AQ. (b) When a multiband 
system with a disconnected Fermi surface has an inter- 
pocket nesting (i.e., strong inter-pocket pair scattering 
and weak intra-pocket one) the gap function has the 
same sign (s- wave symmetry) within each pocket, and 
the nodal lines can happily run in between the pockets. 

Here a natural question is: to prepare multiband sys- 
tems systematically, we can consider lattices comprising 
some units, or "plaquette lattices". This is indeed con- 
sidered by Kimura et al. || who conceived an idea of 
actually fabricating the structure from arrays of quatum 
dots. One example is the inset of Fig. 1, where the unit 
is a square lattice of squres, where the intra-plaquette 
transfer is stronger than the inter-plaquette one. Now, 
an important question is: how the structure of the super- 
conducting gap function having higher T c 's is determined 
and whether the s- wave symmetry on each pocket is com- 
pulsory. 

In this paper, we study the repulsive Hubbard model 
on a square lattice of diamonds, where a diamond is a 
four-site unit rotated by 45 degrees (Fig. 1). One moti- 
vation is that this structure is reminiscent of the in-plane 
lattice of B and C in CaE>2C2 p^ |, so should not be too 
unrealistic, although we are not claiming to consider this 
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particular material. We have solved Eliashberg's equa- 
tion |l3| for the Hubbard model with the FLEX and ob- 
tained a T c that is considerably higher than that of the 
square lattice. More importantly, however, the present 
model, despite of its disconnected Fermi surface, has sign 
changes in the gap function within each pocket, so a dif- 
ferent mechanism should be at work. We can trace back 
its reason in real space that the singlet pairing here re- 
sults from a hybridization of two molecular orbitals in 
each plaquette. 
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FIG. 1. The lattice considered in this paper. The inset 
shows the plaquette lattice in Ref. |J. 

So a message of this work is that, in order to search for 
higher T c 's, we should look not only at the Fermi surface 
but also at the atomic structure in real space. 

We have used the four-band version of FLEX, §,[o|,@ 
where the Green's function G, the spin susceptibility x? 
the self-energy E, and the superconducting gap function 
(j) are 4x4 matrices, e.g., Gi m (k,ie n ), where /,m spec- 
ify the four sites in a unit cell. We can go from the site 
indeces over to band indeces with a unitary transforma- 
tion. For the spin susceptibility, we concentrate on its 
largest eigenvalue denoted as x- From Green's function 
and the spin susceptibility we obtain the superconduct- 
ing gap function (with spin-singlet pairing assumed) and 
T c by solving the linearized Eliashberg's equation [[L3|. 
In our analysis, we take 32 x 32 /c-point meshes and up 
to 4096 Matsubara frequencies, where the numerical re- 
sults are sufficiently converged. The on-site Hubbard re- 
pulsion is taken to be a typical strong-correlation value, 
U/W = 7/6, where W is the band width. The band fill- 
ing (=number of electrons/number of sites) is taken to be 
a value n = 0.85, which is close, but not too close, to the 
half-filling where the Mott transition and antiferromag- 
netic order are expected. In the below we have confirmed 
that when the eigevalue of Eliashberg's equation becomes 
unity prior to the divergence of the spin susceptibility. 

The FLEX result for the spin susceptibility, Green's 
function, and the superconducting gap function are 
shown in Fig. 2 for t' = 0.6 (t = 1.2) with T = 0.06. 
Ridges in the Green's function delineate the Fermi sur- 
face. Fermi surface comprises two bands (second and 
third from the bottom), which we call band A and band 



B. We can see that each band has a squre-like pocket for 
the Fermi surface, which is reminiscent of the Fermi sur- 
face of Sr2Ru04, with a good reason as we shall see. On 
the other hand we notice that the spin susceptibility has 
broad peaks around (7r,7r). Quite unexpectedly, the gap 
function has a strange structure in such as way that its 
amplitude on the Fermi surface is peaked at the corners 
rather than along the edges. 




FIG. 2. Spin susceptibility (x), the absolute value of the 
Green's function squared (|G| 2 ), and the superconducting gap 
function <I> against k x , k y for the lowest Matsubara frequency 
for U = 7, n = 0.85, t 1.2, t' 0.6, and T 0.06. 

Figure 3 shows T c as a function of t. |l5| A plaquette 
lattice is characterized by the intra-plaquette transfer t 
and inter-plaquette one t' . To facilitate comparison, we 
have fixed the single-particle bandwidth W = 2(2|t| + |£ / |) 
to 6, which is the value of W when t = t f = 1. Interest- 
ingly, T c is larger for smaller which is in accord with 
our expectation that the case of plaquette lattice (tetram- 
irization in the present case) is favorable. T c saturates to 
a value T c = 0.06 = 0.01W\ which is two times greater 
than that (< 0.03t) for the square lattice. 

The enhancement of T c in smaller \t'\ can be under- 
stood by the nesting between the disconnected Fermi sur- 
face (Fig. 4). At half filling (n = 1) the Fermi surface is 
perfectly nested regardless of the value of \t'\. When less 
than half-filled, the nesting in the band depicted by Q in 
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Fig. 3 degrades but not so much for smaller \t'\. Because 
the effective attraction here is mediated by the antiferro- 
magnetic spin fluctuation (with a large |Q|), T c is larger 
for better nesting, i.e., for smaller \t'\ ex the warping of 
the quasi- ID Fermi surface. This ^-dependence of T c can 
be understood also in real space, where the intra-cell spin 
singlet should be more robustly formed for smaller \t'\. 
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FIG. 3. T c plotted as a function of t for U = 7 and 
n = 0.85 in units where the band width 2(2t + t') = 6. 
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FIG. 4. Fermi surfaces of the lattice depicted in Fig.l. 
The solid (dashed) curves represent positive (negative) values 
of the gap function on the Fermi surface. The arrows repre- 
sent typical momentum transfer in pair scattering processes, 
where Q corresponds to the wave vector at which the spin sus- 
ceptibility is peaked. The solid (dashed) curves in the inset 
represents the positive (negative) values of the gap function 
in a plaquette lattice in Ref . M . 



We can explain the structures of the strange gap func- 
tion in terms of molecular orbitals (Fig. 5a) and the 
pairing (Fig. 5b) in real space. As mentioned, the Fermi 
surface reminds us of the (a and f3) bands in SrRuO, 
which consists of two sets of quasi- ID array of orbitals. 
In the present case the quasi- ID band structure arises 
from the molecular orbitals in the plaquette. The molec- 



ular orbitals that are relevant to bands A and B have p x - 
and p y - symmetries (Fig. 5a), so that we do have two 
sets of ID arrays of orbitals. Due to the hybridization of 
the two orbtials, the two sets of quasi- ID Fermi surface 
(two, orthogonal sets of parallel lines) anticross, and we 
end up with two square-shaped pockets. 

Now the interesting point is, the Cooper pair, aris- 
ing from antiferromagnetic fluctuations, is formed across 
the adjacent sites (ellipses in Fig. 5b), that is across 
the two molecular orbitals. This implies that the gap 
function on the Fermi surface is large only around the 
positions at which the two (p x and p y ) bands hybridize, 
i.e., around the anticrossing points. The pairings de- 
picted in Fig. 5b as solid ellipses and dashed ellipses 
have opposite signs, so that the symmetry is d xy . We 
have thus a curious case of an inter-band pairing (a pair- 
ing formed by inter-band pair scatterings). We can in 
fact confirm that the pairing is interband and d xy by 
looking at the sign of the gap function. The ellipses 
is Fig. 5 are represented by (Ca(b)^Ca(b)i)i where 
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^4(2)a - CW, where C ia an- 
nihilates an electron with spin a at site i as numbered 
in Fig. 5. When the pairing is formed across adjacent 
sites we can show (Ca^Cai) = — (Cb^Cbi) (opposite 
signs across the bands) with an odd parity (d xy ) against 
k x 4 ^ k x or k y t ^ k y . 

This explains why the gap function changes sign within 
each pocket. Hence the s- wave symmetry on each pocket 
is not a necessary condition for a high T c . Even when 
there are intra-pocket sign changes, the disconnectivity 
of the Fermi surface plays an important role, since the 
inter-pocket sign change associated with Q does exploit 
this in avoiding intersecting the Fermi surface. 

One might then wonder why the intra-pocket sign 
changes do not reduce T c , since the pair-scattering pro- 
cesses depicted by Q' in Fig. 4 connecting the /c-points 
across which the gap function has the same sign would 
normally reduce T c . In the present case, however, we 
have an ingeneous situation where the gap function is 
strongly peaked at the position (corners in this example) 
where hybridization occurs, so that the Q' processes, be- 
ing off the peak, have little effect. 

How is the present model compared with the lattice 
considered in Ref. ||, where the plaquette is a four- 
membered ring in either case. Here the consideration 
in real space comes in handy. The latter case is a 2D sys- 
tem rather than quasi-lD, where the pockets arise due 
to a band folding due to the tetramerization,. The pair- 
ing in real space is formed within each molecular orbital 
(Fig. 6) unlike the present case, so the superconducting 
gap function has an s-wave symmetry on each Fermi sur- 
face pocket (but opposite signs across the two pockets 
with the nesting Q w due to the band folding (inset of 
Fig. 4). Despite of this distinction, the resulting T c 's are 
similar between the two cases. 
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Another interesting comparison is with S^RuC^ 
[|l6],[L7|, which also suggests an importance of the real 
space structure. Although the Fermi surface is quasi- ID 
in Sr2Ru04 as well, the oxide has strong pieces of evi- 
dence for triplet pairing, and even when one considers the 
spin-singlet pairing, T c (for the a, f3 bands for this ma- 
terial) estimated with the FLEX is very small |l7| ]. Due 
to the quasi-one dimensionality the spin susceptibility 
x(QiW) has linear ridges in k-space. As a result, the pair 
scatterings have large contributions all over the ridges, 
so that the (extended) s-wave gap function involves un- 
favorable pair scatterings across which the gap has the 
same sign, resulting in the reduced T c . This makes the 
spin-triplet, p-wave pairing more favorable, for which the 
resulting gap function has large amplitudes all over the 
Fermi surface, except for dips at the corners . By con- 
trast, the present case has a built-in antiferromagnetic 
structure within the unit cell in real space, so that the 
spin susceptibility has a well-defined peak around (7r,7r). 
This makes a specific Q to be relevant in the pair scat- 
tering, which makes the gap function peaked at specific 
points, which in turn gives rise to an enhancement in T c . 

To summarize, we should question not only the struc- 
ture of the Fermi surface but also the structure of the 
molecular orbitals and the superconducting pairings in 
real space in understanding the superconductivity from 
the electron repulsion. One interesting tendency is that 
high T c 's are found in plaquette systems where unit cells 
having multiple sites are connected with a relatively small 
inter-plaquette hopping [p|-p^0|, which also suggests the 
importance of the real space picture. In fact, the rela- 
tion between real-space and the momentum-space pic- 
tures has been discussed for systems consisting of dimers 
recently. |10| It is an interesting future problem to see 
how the real-space and the momentum-space pictures are 
related with each other for wider class of systems in re- 
alizing high T c 's. 
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FIG. 6. (a) Intra-plaquette molecular orbitals with p x (al) 
and p y (a,2) symmetries, and (b) Cooper pairs for the lattice 
considered in Ref. [Hi. 
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FIG. 5. (a) Intra-plaquette molecular orbitals with p x 
(al) and p y (a2) symmetries, (b) Cooper pairs with solid and 
dashed ellipses indicating pairing amplitudes with opposite 
signs. 
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